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Abstract 

New classes of exact solutions of the three-dimensional unsteady Navier-Stokes 
equations containing arbitrary functions and parameters are described. Various 
periodic and other solutions, which are expressed through elementary functions 
are obtained. The general physical interpretation and classification of solutions is 
given. 
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1 Class of motions of the viscous incompressible fluid 
under consideration 

Self-similar, invariant, partially invariant, and certain other exact solutions of the 
Navier-Stokes equations including those with generalized separation o f variables were 
considered, for example, in [1-15]. Below, the term "exact solutions" is used according 
to the definition given in [14, p. 10]. 

Three-dimensional unsteady motions of a viscous incompressible fluid are de- 
scribed by the Navier-Stokes and continuity equations: 

~aT + Vl ~x~ + V *~x~ + ^3^— = -zr + v \ -ET + TT" + "TT" > W 

dt ox dy dz p dx \ ox 2 - dy 2 - dz z J 

dv 2 { av 2 | av 2 | v dv 2 _ idP [ Jd 2 v 2 | d 2 v 2 | d 2 v 2 \ 

dt 1 dx 2 dy 3 dz p dy \ dx 2 dy 2 dz 2 J ' 

^Yl + v ^ + V — + V — - --— + v (^Yl + ^Xl + (3) 
dt 1 dx 2 dy 3 dz p dz \ dx 2 dy 2 dz 2 J ' 



dVx dV 2 dV 3 
dx dy dz 



= 0. (4) 
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Here x, y, and z are the Cartesian coordinates, t is time; V\, V%, and V3 are the fuid- 
velocity components; P is pressure; and p and v are the fluid density and kinematic 
viscosity, respectively. When writing Eqs. (l)-(4), it was assumed that the mass forces 
are potential and included in the pressure. 

We consider the flow of a viscous incompressible fluid when the fluid-velocity 
vector on the z axis is directed along this axis. Near the z axis, the transverse velocity 
components are small and can be expanded in a Taylor series in terms of the transverse 
x and y coordinates. If we restrict ourselves to the main terms of the expansion in x 
and y for the velocity components, it is possible to obtain the following representation 
for the desired values after the corresponding analysis: 

/ 1 OF \ / 1 dF \ 

Vl = X (--—+w)+yv, V 2 =xu-y(-—+w), V 3 = F, 

p 1 2 1 o 2 1 ^2 dF f dF , ^ 

— =Pa--ax --py -jxy--F +v- / -^-dz, 

p 2 2 2 dz J at 

where pg, a, (3, and 7 are arbitrary functions of time t setting the transverse pressure 
distribution; F, u, v, and w are unknown functions dependent on the coordinate z and t. 
The substitution of Eqs. (5) into Navier-Stokes Eq. (3) and continuity Eq. (4) results in 
identities, and Eqs. (1) and (2) become A n x + B n y = (n = 1, 2), where A n and B n 
represent certain differential expressions dependent on the variables z and t. The split 
in the variables x and y results in four equations A n = and B n =0 (n = 1, 2), 
which can be transformed to the following form: 

d 2 F d 2 F l/dF\* , n , d 3 F . 2n 

+ ^-5-5- (-5-) = + +"tt + 2{uv + w 2 ), (6) 



dtdz dz 2 2 V dz / 
du „ du dF 
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at oz oz oz 1 * 
dv dv dF d 2 v 

^- +F— -v— = l + vjpt , (8) 
Ot Oz oz OZ z 

dw ^ pdw dF a — (3 ^ d 2 w 

dt dz dz 2 dz 2 

It is important to emphasize that solution (5) precisely satisfies Eqs. (l)-(4) for the 
viscous fluid motion by virtue of Eqs. (6)-(9). 

For 7 = 0, the structure of exact solution (5) and system (6)-(9) was obtained in 
[12] from other reasons by the investigation of the class of partially invariant solutions 
(the case ofa = /3 = 7 = was considered in [7]). In [12], the group classifica- 
tion of system (6)-(9) was carried out for 7 = 0, which resulted in singling out two 
types of time dependences for the determining functions: (i) a and /3 are constant, and 
(ii) a and ft are proportional to t~ 2 (the exact solutions of the Navier-Stokes equations 
with a reasonably simple structure correspond to these dependences). 

In this work, we obtained new classes of exact solutions of system (6)-(9), when 
the determining functions contain a functional arbitrariness. The basic idea of the 
following analysis is that we can obtain a single isolated equation for the longitudinal 
velocity component V3 = F from system (6)-(9). 
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2 Reduction of system (6)-(9) to a single equation 



We consider first the special class of exact solutions described by a single equation. 
In Eqs. (6)-(9), we put 

dF dF dF 

u = m— + A, v = n— + B, w = k— + C, (10) 

oz oz oz 

where m, n, k, A, B, and C are the desired functions of time t. We require that 
four Eqs. (6)-(9) coincide after the substitution of Eq. (10) in them. As a result, for 
determining the desired functions, we obtain the nonlinear system consisting of one 
algebraic equation and six ordinary differential equations: 

9 1 

mn + k = -, (11) 

A_m^ = B_j± = C-k[ = +Bm + 

m n k 

a - ir,' 

-a- f3 + 2AB + 2C z . (13) 



m 



2k 



This system contains seven equations for nine functions — six functions m, n, k, A, B, 
and C from Eq. (10) and three functions a, f3, and 7 from Eqs. (6)-(9) (in this case, they 
are also treated as desired). It is possible to show that the last equation in Eq. (12) is 
the consequence of three other equations (1 1), (12). Therefore, three desired functions 
in system (1 1), (13) can, in general, be chosen arbitrarily. 

Taking into account Eqs. (10)— (13), we reduce system (6)-(9) to a single equation 

+F— - (—) 2 - v— +q— +p (14) 

dtdz dz 2 V dz J dz 3 dz ' 

the functions p = p(t) and q = q(t) are defined by the relations 

1-A' A-m! 

P= , 9= • (15) 

m m 

The general property of Eq. (14). Suppose Fq(z, t) is a solution of this equation. 
Then the function 

F = F (z + i>(t),t)-^ t (t), (16) 

where ip(t) is an arbitrary function, is also a solution of Eq. (14). 

For constructing solutions of system (11)-(13), it is necessary to distinguish two 
cases. 



1°. Case ofm = n. In this case, the general solution of system (1 1)— (13) can be 
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represented as 

m = n = g sm ^) & = h cos V) 

A = B = ^(q sin 93 + <^ cos ip), C = i (g cos tp — <p' t sin <p) , 

" = j? 2 + i(^) 2 - fKl - cos^) + CJ, (17) 

13 = \q 2 + W t ?-\p(l + cos p)-C' t , 

7 = ipsin^ + AJ, 

where p — p(t), q = q(t), and <p = ip(t) are arbitrary functions. For convenience, 
the free functions p and q in Eq. (17) are chosen so that system (6)-(9) is reduced to a 
single equation (14) with the same functions p = p(t) and q = q(t) as a result of the 
transformation (10), (17). 

Thus, this important statement is proved. An arbitrary solution of Eq. (14) for 
arbitrary functions p = p(t) and q = q(t) generates an exact solution of the Navier- 
Stokes equations (l)-(4). This solution is described by the function F = F(z, t) and 
Eqs. (5), (10), (17). 



2°. Case of m ^ n. In this case, the general solution of system (11)-(13) can 
be obtained as follows. The functions m = m(t), k = k(t), and q = q(t) are set 
arbitrarily under the condition m 2 + k 2 ^ 1 /4. The remaining functions included in 
system (1 1)-(13) and Eq. (14) are calculated sequentially from the formulas 



n 



(18) 



4m 

A = mq + m! t , B = nq + n' t , C = kq + k' t , 
P = A ' t ~ B ' t , a = AB + C 2 + C{-l P (l-2k), 
0= AB + C 2 -C' t - -p(l + 2k), j=pm + A' t . 

In this case, the coefficient p = p{t) in Eq. (14) is determined through the functions 
m = m{t), k = k(t), and q = q(t) and their derivatives (instead of being set arbitrarily 
as in the case of m = n). An attempt to set p = p(t) directly instead of the function m 
(or k) results in a nonlinear ordinary differential equation of the second order for the 
function m (or k) with an arbitrary function q = q(t). 

We consider how to choose the function q = q(t) so that the identity p = is 
satisfied. From the expression for p in Eq. (18), we have A = B + sq, where so is an 
arbitrary constant. From here, taking into account Eqs. (18) for A, B, and n, we find 
the function 

4(m z + Ar) — 1 m 
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3 Exact solutions of Eq. (14) for various p = pit) and 
q = q(t) 



1°. Functional separable solution: 

F = -a' t (t) + b(t)[z + a(t)} - z ^ , q = -4b, P = b' t + 3b 2 , 
where a = a(t) and b = b(t) are arbitrary functions. 

2°. Periodic solutions in the form of the product of functions of different argu- 
ments: 



F = a(t) s'm(az + B), a(t) = C cxp[-j/cr 2 t + J q(t) dt 
p = —a 2 a 2 (t), q = q(t) is an arbitrary function, 



(20) 



where B, C, and a are arbitrary constants. Putting in Eq. (20) that q(t) = va 1 + ip' t (t), 
where (p(i) is a periodic function, we obtain a solution periodic in both arguments z 
and t. 

Example. Consider the stationary case. In Eqs. (17), (20), we put 

i-p = 0, a = ^— — -, q = va 2 = 2a\, p=—a 2 a 2 , a = {2a\jvfl 2 . 

a 

As a result, using Eqs. (5) and (10), we obtain the solution 

Vi = a%x, V2 = \(ai + a 2 ) cos(ctz) — ai]y, Vs = sin(a-z), 

a 

which describes the three-dimensional flow of a fluid layer between two flat elastic 
films (the film position depends on the values of z = and z = 2-k/<j), the surfaces of 
which are stretched according to the law V\ = a\x and V% = a-iy- 



3°. Generalized separable solutions exponential in z: 

F = a(t)e^ IJZ + b(t), p = 0, q= — - ab - a 2 v, (21) 

a 

where a = a(t) and b = b(t) are arbitrary functions. Choosing a(t) and b(t) to be 
periodic functions, we obtain a solution periodic in time. 

Formulas (20) and (21) together with Eqs. (5), (10), (17) define two classes of 
solutions of the Navier-Stokes equations dependent on several arbitrary functions. 
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4°. Solution (21) can be represented in the form 



-az + a 2 vt + (q + ab)dt +b(t), p = 0, (22) 



F = clq exp 

where b = b(t) and q = q(t) are arbitrary functions and ao is an arbitrary constant. 
Formula (22) defines a new class of exact solutions of the Navier-Stokes equations 
with the help of Eqs. (5), (18), forp = and Eq. (19). 

5°. Exact solution in the form of the product of functions of different arguments: 

F = a(t)(C 1 e' JZ +C 2 e- az ), p = 4CiC 2 <7 2 a 2 (t), q = ^ - o 2 v, 

a 

where a = a(t) is an arbitrary function, C\, C2, and a are arbitrary constants. 

6°. Monotonic traveling-wave solution: 

F = -Qua tanh[cr(z - Xb) + B] + X, p = 0, q= &vo 2 . 

7°. Unbounded periodic traveling-wave solution: 

F = &ua tan[cr(z -Xt)+B]+X, p = 0, q= -%vo 2 . 

8°. Functional separable solution: 

where the functions a = a(t), b = b(t), c = c(t), and A - X(t) satisfy the system of 
ordinary differential equations 

X' t = — cA, a' t = (v\ 2 + q + 2c + bX)a, c' t = c 2 + qc + p. 

Here, three of the six functions a(t), b(i), c(t), X(t), p(t), and q(t) can be set arbitrarily. 

9°. Functional separable solution: 

F = u(t)z + ^sm[9(t)z + a], (23) 

where a is an arbitrary constant, and the functions u> = uj(t), £ = £(t), and 6 = 6(t) 
are described by the system of ordinary differential equations 

tf t = -uB, u' t = oj 2 + q{t)uj + p(t) + f , ^ = [2w - v6 2 + q(t)]£. (24) 

In this system, it is possible to treat the functions 0(t) and as arbitrary, whereas 
the functions w(t), p(t), and q{t) are elementarily determined (without quadratures). 
A periodic solution (23) corresponds to periodic functions 6{t) and 
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10°. Functional separable solution: 



C ie 8 ^ + C 2 e-°^ 



(25) 



where C% and C2 are arbitrary constants, and the functions ui = ui(t), £ = £(t), and 
6 = 6(t) are described by the system of ordinary differential equations 



-lo9, w't = u 2 + q{t)u+p{t) - 4CiC 2 f , £ = [2w + ^ 2 + g(t)]£. (26) 



Remark. See also [8, 13, 15] for exact solutions of Eq. (14) with q = 0. 



4 Reduction of system (6)-(9) to two equations 

We describe two cases of reducing system (6)-(9) to a single isolated nonlinear equa- 
tion for the longitudinal velocity F and a second equation for determining a new aux- 
iliary function. 



1°. First case. By letting 

u = a 2 G, v = ~b 2 G, w = abG, a = 0, 7 = 0, (27) 

where a and b are arbitrary constants, we reduce system (6)-(9) to an isolated equation 
for the longitudinal velocity F and an additional equation for the function G = G(z, t): 

d 2 F d 2 F 1/&FV „ d 3 F 

— + F— - G— - u— (29) 

dt dz dz dz 2 



2°. Second case. In Eqs. (6)-(9), let 



1 • ( dF \ 2r, 

u = - sin 03^— + qj + a B, 

1 ( dF \ 

w = - cos 09 — h q + ab<d, 

2 V / 

1 2 1 /, x 1 / 

""I 9 _ 2 P ^ ° S ^ + 2 9tCOS</: '' 

/3 = \q 2 - + cos <p) - -q' t cos ip, 

1 1 / • 

7 = -p sin 95+ -g t sino3, 



(30) 
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where p = p(t) and q = q(t) are arbitrary functions, and a and b are arbitrary con- 
stants, 9 = Q{z, t) is an unknown function, and tp is the constant determined from the 
transcendental equation 

(a 2 - b 2 ) sinip + 2abcosip = 0. (31) 

As a result, system (6)-(9) is reduced to two equations: 

d 2 F +p d 2 F f dF \ 2 _ d 3 F + dF + ^ 
dtdz dz 2 \ dz J dz 3 dz 

oe oe bf o 2 e 

at oz oz oz z 

The nonlinear equation (32) for F coincides with Eq. (14) and can be treated indepen- 
dently (some of its exact solutions were described previously), and Eq. (33) is linear 
with respect to the desired function Q. 

For stationary solutions of Eqs. (32) and (28) (for constant p, q, and a), the nonsta- 
tionary equations (33) and (29) are linear separable equations, whose solutions can be 
obtained using the Laplace transform in time. 

Equation (32) (and Eq. (28)) admits an obvious degenerate solution F = a(t)z + 
b(t); in this case, the corresponding Eq. (33) (and Eq. (29)) can be reduced to the linear 
heat equation. 

System (28), (29) for an arbitrary function a = a(t) has a solution in the form 

F = az 2 + b(t)z + -^\b 2 {t) - 2b' t (t) - 4a(t)], 
w = A{t)z 2 + B{t)z + C(t), 

where a is an arbitrary constant (a ^ 0), b(t) is an arbitrary function, and the func- 
tions A = A(t), B = B(t), and C = C(t) are described by the system of ordinary 
differential equations, which is not presented here. 



5 Interpretation and classification of the flows under 
consideration 

Arbitrary fluid flows having two symmetry planes admit a representation of the type 
of Eq. (5) in the vicinity of the line of intersection of these planes (in the adopted 
notation, the planes intersect in the z axis). Such flows include the axisymmetric flows, 
combinations of axisymmetric flows with rotation around of the z axis (in particular, 
the von Karman type flows), plane flows symmetric with respect to a straight line, 
flows in rectilinear impenetrable and porous pipes with elliptic and rectangular cross 
sections, fluid jets flowing from orifices of elliptic and rectangular shapes, etc. (see also 
[10,11]). 

It is convenient to treat the axial flows described by Eqs. (5) as a nonlinear super- 
position of a translatory (nonuniform) flow along the z axis and a linear shear flow of 



a special type. In the vicinity of the point z = zq lying on the axis, the components of 
fluid velocity taking into account Eq. (5) can be represented as 

Vk — F$k3 + Gk m X rn ; 
G lx = -\F z + w, G 12 = v, G 21 =u, G 22 = -\F z -w, 

(34) 

G13 = G23 = G31 = G32 = 0, G33 = F z \ 
Xi = x, X 2 = y, X 3 = z-z . 

Here k, m = 1, 2, 3; the Gkm are shear matrix components; the summation is as- 
sumed over the repeated subscript m; 6km is the Kronecker delta; and F z is the partial 
derivative with respect to z. All values in Eq. (34) are taken for z = zq. The vanish- 
ing of the sum Gn + G22 + G33 = of diagonal elements is a consequence of fluid 
incompressibility. 

An arbitrary matrix 1 1 Gkm | can be represented in the form of the sum of a symmet- 
ric and an asymmetric matrix 

||Gfc m || = H-Efemll + ||Ofc rn ||, 

1 1 (35) 

Fkm ^mfc 2 {Gkm G 7n k)i ^km ^mfc 2 {Gkm G m ^J. 

In turn, the symmetric matrix ||-Efc m || (in this case, it can be identified with the 
strain-rate tensor) can be reduced to a diagonal form with diagonal elements E%, E 2 , 
and E3, which are roots of the cubic equation dct \\Ek m — A<5fc m || = for A, by 
appropriately rotating the system of coordinates. 

For this flow (34), the diagonal elements determining the intensity of the tension 
(compression) motion along the respective axes are calculated from the formulas 



£1,2 = ~^F Z ± - y/4w 2 + {u + v) 2 , E 3 = F Z . (36) 

The splitting of the shear coefficient matrix 1 1 Gkm 1 1 into symmetric and asymmetric 
parts (35) corresponds to the representation of the velocity field of the linear shear 
flow of the fluid as a superposition of a linear deformational flow with the tension 
coefficients E\,E 2 , and £3 along the principal axes and rotations of the fluid as a solid 
body with the angular velocity uj = (^32, SI13, ^21)- 

For this flow (34), we have ^32 = fiia = and the fluid rotates around the z axis 
with the angular velocity 

n 21 = ^(u-v). (37) 

It is easy to show that Eqs. (36) and (37) remain valid for an arbitrary point 
(xq, yo, Zq) of the flow (5) under consideration. 

The analysis of Eqs. (36), (37) enables us to single out certain characteristic types 
of flows indicated in the classification table. 
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Table 1 : Classification of axial flows described by Eqs. (5) 



Type of flow 


Desired functions 


Functions included in pressure 


Axisymmetric 


u = v = w = 


a = /3, 7 = 


Combination of axisymmetric flow 
and rotation around of the z axis 


ui = 0, v = —u 


a = /3, 7 = 


Pure deformational (without rotation) 


u = V 


a, 0, 7 are arbitrary functions 


General axial 


u 5^ v 


a, /3, 7 are arbitrary functions 



6 Some generalizations 

Let Vi(x, y, z, t), V^x, y, z, t), Vs(x, y, z, t), and P(x, y, z, t) be a certain solution of 
Navier-Stokes equations (l)-(4). Then the set of functions 

Vi = Vi(x-x , y - y , z - z ,t) +x' , 
V 2 =V 2 (x- x ,y- y ,z-z ,t)+y' , 
V3 = V 3 (x -x ,y- y , z - z , t) + z' , 
P = P(x -x ,y- yo,z- z ,t) - p(x'qX + y'^y + z'qz), 

where xq = xo(t), yo = yo(t), and zq = zo(t) are arbitrary functions (primes de- 
note the derivatives with respect to t), also gives the solution of Eqs. (l)-(4) [4, 15]. 
The combination of Eqs. (5) and (38) for zq = determines an exact solution of the 
Navier-Stokes equations, which can be treated as the generalized axial flow with the z 
axis moving along the plane x, y, according to the law x = xo(t), y = yo(t). The indi- 
cated solution can be used for the mathematical simulation of destructive atmospheric 
phenomena such as waterspouts and tornados. 
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